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Abstract 



^ I The two-frequency sine-Gordon model is examined. The focus is mainly on the 
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nonperturbative (truncated conformal space approach) framework, and give particu- 
lar attention to a phase transition conjectured earlier by Delfino and Mussardo. We 
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Ising universality class. Furthermore, we check the UV-IR operator correspondence 
and conjecture the phase diagram of the theory. 

PACS codes: 64.60.Fr, ll.lO.Kk 

Keywords: nonintegrable field theory, two-frequency/double sine-Gordon model, form fac- 
tors, truncated conformal space approach, finite size effects, phase transition 



*bajnok@poe. elte.hu 
^palla@ludens. elte.hu 
■''takacs@mth.kcl.ac.uk 
§ wf eri @ af avant .elte.hu 



1 



1 Introduction 



The sine-Gordon model has attracted a great deal of interest over the decades, for the par- 
ticular reason that while being an integrable field theory, it can be used as a toy model for 
many nonperturbative quantum field theory phenomena. It also has very interesting appli- 
cations in diverse areas of physics, ranging from statistical mechanics of one-dimensional 
quantum spin chains to nonlinear optics (for a - non-exhaustive, but representative - list 
with references see the introduction of ||l|). 

In the present paper we study a nonintegrable extension of sine-Gordon theory in which 
the scalar potential consists of two cosine terms with different frequencies, called the two- 
frequency or double sine-Gordon (DSG) model. It is suggested in that this model can 
be used as a more refined approximation to some of the physical situations (e.g. wave 
propagation in a nonlinear medium) where the ordinary sine-Gordon model is applicable. 
The model is also interesting theoretically, since due to its nonintegrability it is expected 
to have more general behaviour than usual sine-Gordon theory and therefore it can be seen 
as a more realistic toy model for nonperturbative quantum field theory. Applications to 
the study of massive Schwinger model (two-dimensional quantum electrodynamics) and a 
generalized Ashkin-Teller model (a quantum spin system) are discussed in ||l] and another 
potentially interesting application to the one-dimensional Hubbard model is examined in 
(together with the generalized Ashkin-Teller model mentioned above). A further po- 
tentially interesting application of the two- (and multi-) frequency sine-Gordon model is for 
ultra-short optical pulses propagating in resonant degenerate medium 0. 

We show that for a rational ratio of the frequencies this model is still amenable to a 
nonperturbative treatment using the truncated conformal space approach (TCSA), and 
that a great deal of interesting information can be extracted, in particular the existence 
and the characteristics of a phase transition which was predicted in on the basis of 
simple classical (mean field) arguments. 

The outline of the paper is as follows. In Section || we introduce the model and discuss 
some of its basic features. Then in Section ^ we set up the TCSA framework for the 
model. Section || discusses results obtained from the so-called form factor perturbation 
theory (FFPT), while in Section |5| we compare the FFPT results to the numerical data 
resulting from TCSA in order to test the reliability of the latter method. Section ^ is 
devoted to a preliminary discussion of phase transition in the case when the ratio of the 
two frequencies is 1/2 (denoted by DSG2) and contains a general analysis of signatures 
of second and first order phase transitions in finite volume in a separate subsection. In 
Section we apply these general considerations to construct the phase diagram of the model 
DSG2 and then we make our conclusions in Section |^. The paper contains an Appendix 
which collects the formulas used for computing the first order corrections in form factor 
perturbation theory. 

2 The two-frequency sine-Gordon model 
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2.1 Definition of the model 

The two-frequency sine-Gordon (or double sine-Gordon, DSG) model is defined by the 
action 

Absg = j dt J dx Q9^$5'^$ + /icos/3$ + Acos(a<l> + 5)^ , (2.1) 

where /3, a are the two frequencies, /i, A are dimensionful couplings (classically of dimension 
mass^) and 5 is a relative phase. Setting yU or A to zero the DSG model reduces to an 
ordinary sine-Gordon model. 

At the quantum level this theory can be considered as a double perturbation of its c = 1 
UV limiting conformal field theory. The action can be written in the following form 

-^DSG = Ac=l — Apert , (2.2) 

where 

Ac=i = j dt j dx^a^$a'^$ , (2.3) 

^pert = -\jdtjdx {fiVf, + fiV-p + Ae*V„ + Ae-^V_„) , (2.4) 

and denotes the exponential field (vertex operator) 

= e^"* , (2.5) 

which in the UV limit corresponds to a {U{1) Kac-Moody) primary field with conformal 
dimensions 

A^ = A^ = f, (2.6) 

where the upper index ± corresponds to the right/left conformal algebra. Correspondingly 
the couplings at the quantum level have the following dimensions 

[/i] = (mass)'-'^^ , (2.7) 
[A] = (mass)'-'"^" . (2.8) 

We also introduce two mass scales and Mg, which are defined to be the masses of the 
sine-Gordon solitons in the models with /i = or A = 0, respectively. It is known ||^ that 
these masses are related to the couplings A, fi via 

/i = K^M^ , 

A = K^Ml-^^\ (2.9) 

where 



2r(A^) I v^r 

7rr(l - A^) 2r ( 




(2.10) 



2-2A 
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2.2 General properties of the model 

The DSG model was examined in by Delfino and Mussardo, where some general obser- 
vations were made. Here we recall these facts without eventually going into any details. 

The first observation is a simple consequence of the formulation of the model as a 
perturbed c = 1 CFT. Namely, in order for the two perturbing operators to be relevant, 
one has to restrict 

f3^<8TT , a^<87i. (2.11) 

Furthermore, we assume that the theory is renormalizable in the strict sense, i.e. all 
divergences can be absorbed into redefining the coupling constants and A (otherwise 
one needs to add new (periodic) counterterms to the action and the theory is no more a 
two-frequency sine-Gordon model but rather some multi-frequency version). This means 
that we have to restrict 

a(3<47T. (2.12) 

In the sequel we always assume that the above two conditions hold. 

Another nontrivial property of the model is that it is a nonintegrable field theoryQ if 
a ^ [3 and both couplings A, /i are nonzero. In addition, the model exhibits drastically 
different behaviour depending on whether the ratio of the two frequencies is rational or 
irrational. 



(2.13) 



2.2.1 Rational case 
If 

a m 

with m and n two relative prime positive integers, then the potential — /i cos A cos (a$ + 6) 
is periodic with period 

2™ 27rm 

— = — • <"^' 

It can then be shown that the fundamental range of the 6 parameter is 

\5\ < I (2.15) 

(i.e. any model with other values of the parameter 6 can be redefined into one with 
a 5 lying in the fundamental range by shifting the field with some integer multiple of 
^). The periodicity of the potential implies that the theory has an infinitely degenerate 
ground state and supports topologically charged excitations. The fundamental topological 
excitation degenerates in the A = limit to an n-soliton state of the corresponding sine- 
Gordon theory and similarly for the /i = limit it is an m-soliton state. For general values 
of the couplings, the solitons are in a sense "confined" into these composite objects. 

^One can argue for the nonintegrability e.g. noting that the so-called "counting argument" fails to 
provide any conserved higher spin quantities \{ a ^ (3 and both A and /i are nonzero. 
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2.2.2 Irrational case 



If, on the other hand, 

^^Q, (2.16) 

then the potential is not periodic. There are no topologically charged excitations and 
the solitons are completely confined. Delfino and Mussardo ^ argue that in this case 
the ground state is unique. Furthermore, from (|2.15|) we can infer that the parameter 6 



presumably plays no role in the dynamics of the model. 

The irrational case is extremely complicated, even as a quantum mechanical problem 
(i.e. neglecting the x dependence). In the next section we shall see that it is impossible to 
apply the well-known TCSA method in this case and in lack of any alternative nonpertur- 
bative method in the present paper we restrict our attention to the rational case (although 
the perturbative formulas presented in Section ^ apply in the irrational case as well). 



3 Truncated Conformal Space Approach 

Here we only review the specifics of TCSA as applied to the DSG model. This nonpertur- 
bative (numerical) method was invented by Yurov and Zamolodchikov in ||5| and extended 
to c = 1 theories in to which we refer the interested reader for more details. 



3.1 The UV Hilbert space 

The first important issue is the operator content of the UV limiting c = 1 CFT. We need 
a theory in which the perturbing operators V±a, y±i3 correspond to local operators. 

In the c = 1 CFT we consider the spectrum of local operators is generated by the 
conformal families of U{1) Kac-Moody primary fields Vr,s with scaling dimensions 

<s = \(^^^\sR^ .r.se'L, (3.1) 
where R is the compactification radius of the c = 1 free boson field x with the action 

^x = ^j ^^^d^xd^x , X-X + 2nR. (3.2) 

Taking into account eqn. (|2.13| ) and the fact that the normalization of $ and x differ by 
\/47r, their identification is essentially unique up to a positive integer k: 



„ , v47rn , v47rm 

R= k-^-— = k- , 3.3 

p a 

y±l3 = V±nfc,0 ) = V±mfc,0- (3.4) 
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The integer k is the folding number, which has been considered in the context of sine- 
Gordon theory in 0] and can be introduced in any scalar field theory with periodic po- 
tential: it corresponds to identifying the field as an angular variable after k periods of the 
potential. We choose k = 1, which means we identify 

$ = $ + —-/ , leZ. (3.5) 

However, it can now be seen that in the limit A = we recover an n-folded version of 
sine-Gordon theory (denoted SG(/3,n) in 0) while in the /i = case we get the theory 
SG(a,m). We remark that by choosing A; = 1 we do not lose any information on the DSG 
model, since the other k > 1 folded versions can be obtained from the k = 1 case by adding 
some twisted sectors as in 0. However, we shall see that the fact that the sine-Gordon 
limits are in general folded models will play a very important role in the dynamics of the 
DSG model. There is also a fermionic version of the theory which has primary fields Vr,s 
in the UV with 

seZ,reZ+^ (3.6) 

but similarly to the folded ones, this version can be constructed by replacing some of the 
topologically charged sectors with some (fermionic) twisted ones. As we are only interested 
in the dynamics of the vacuum sector we shall not discuss this modification here either. 

We also wish to note that as the irrational case can be considered as a limit n, m oo, 
the Hilbert space of the irrational DSG model corresponds to an uncompactified boson in 
the ultraviolet limit, with a continuous spectrum of primary fields 

2 

V.,0 , ^? = y , reR. (3.7) 

The fact that the spectrum is continuous means that TCSA is not applicable since the 
conformal Hilbert space contains infinitely many states even when truncated by imposing 
some upper bound on conformal energy. 



3.2 The TCSA Hamiltonian 

We put the system on a cylinder with finite spatial volume L, i.e. we take the space 
coordinate < x < L and identify x with x + L. Given the identification (|3.5|) this means 
in fact the following quasi-periodic boundary condition for the field 

^(x + L) = <l>(x) + — — g , qeZ. (3.8) 

q is the winding number labelling different sectors of the theory, and corresponds to a 
conserved topological charge in the full perturbed CFT. In the present paper we are only 
interested in the g = sector as it contains all the relevant information for the problems 
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treated in this work. In addition, in the sector containing the vacuum the Lorentz spin of 
all the states is also zero, i.e. (Lq — Lq) |\[') = for any state 

After a mapping onto the conformal plane, the Hamiltonian of the theory can be written 
in the following form 

(3.9) 

The matrix elements of this operator can be calculated explicitly between any two 
states in the UV limiting Hilbert space. Truncating this space to finitely many state by 
introducing an upper conformal energy cutoff the Hamiltonian becomes a finite matrix, 
which can then be diagonalized numerically to get an approximation of the spectrum. 
Thus one can extract the energy E^{L) of any state \l/ as a function of the volume L. 

It is well-known P| that TCSA suffers from UV divergences whenever the scaling di- 
mension A"*" + A~ of any of the perturbing operators is larger than 1. In this case the 
only quantities that are possible to extract are the energies relative to the ground state 
i.e. Ei{L) — Eo{L). Furthermore, in the case where TCSA is UV divergent the truncation 
errors are much larger even for such relative energies, and grow very fast with the volume 
L, which restricts the usefulness of TCSA in this case to small values of L. 



3.3 The interpolating mass scale 

In order to get numerical results, we have to fix the units in which we measure energies 
and distances. For later convenience, we introduce the following 'interpolating' mass scale 

M = M££l±^, (3.10) 

where = 2 — 2A^. We also use the following dimensionless combination of the coupling 
constants 

V= . ^x., • (3.11) 

The advantage of this parametrization is that M interpolates smoothly between the mass 
scales of the two limiting sine-Gordon models, which correspond to ?7 = and 1, respec- 
tively. Introducing the notation 

1 1 1 1 1 1 

D=^- + ^^-^, (3.12) 



we obtain 



Ai= ^ , A = ^' , 3.13 
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and the Hamiltonian can be made dimensionless as follows: 



1 



+ , (3.14) 



2(27r)^'--i V^, 

where we introduced the dimensionless volume parameter 

1 = ML. (3.15) 

We denote the eigenvalue of /^dsg corresponding to the state and as a function of / by 
e^{l) i.e. E^{L) = Me^{ML). 

It is very important that /^dsg interpolates smoothly between the two extremal sine- 
Gordon models as t] varies from to 1. Had we chosen as our units one of the soliton 
mass scales Mp or M^, the dimensionless Hamiltonian would have been singular at the 
other end of the rj interval. We shall discuss that there is a little inconvenience in the 



choice (|3.10|) when we want to compare with perturbation theory around one of the end 
points, since such a calculation is more naturally formulated in terms of the mass scale of 
the corresponding extremal sine-Gordon model rather than the interpolating mass scale. 
However, the convenience of being able to scan the whole range of rj with a single numerical 
program well outweighs this disadvantage. 



4 Form factor perturbation theory 
4.1 Generalities 

To supplement and check TCSA (which is a numerical method) with some analytic results, 
we use the so-caWed form factor perturbation theory (FFPT). The viewpoint of considering 
a nonintegrable field theory as a perturbation of an integrable one was first taken in to 
which we refer the reader for a detailed exposition of the subject. Here we restrict ourselves 
to a brief summary of the necessary formulae. 

Consider an integrable quantum field theory with action perturbed by a spinless 
local field "^{x): 

A = Ao-g j d?x^{x) . (4.1) 

Assume that the spectrum of the unperturbed theory is known, together with the 
following matrix elements of \l/(x) (form factors): 

FZ...a^ (^1, . . . , = o(0|vl/(0)|ai ...an (^„))|r , (4.2) 
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where |0)g is the vacuum state and |ai . . . a„ (^9n))o° denotes an asymptotic n-particle 
in-state of particles oi, . . . , with rapidities . . . , /9ri in the unperturbed theory ho- 
using FFPT one can calculate the spectrum of the theory A perturbatively in the coupling 
constant g. Here we list the results obtained to first order. The vacuum energy density is 
shifted by an amount 

5^vac = ^7o(0|^|0)o. (4.3) 

The mass (squared) matrix M^^ gets a correction 

5Mf^ = 2gFZ (m , 0) 5^^,^, (4.4) 

supposing that the original mass matrix was diagonal and of the form 

Mf, = m%b . (4.5) 

Finally, the scattering amplitude for the four particle process a + b — > c + c? is modified by 

F*- , (in, t9 + in, 0, t9) 
niami) smh 

4.2 Analytic results 

Using the general formulas above and the ones for the form factors and vacuum expectation 
values, listed in Appendix one can derive analytic formulas for the first order corrections 
to the vacuum energy density and the spectrum. 

We identify the unperturbed theory and the perturbing operator in the following way: 

Aq = j j Qa^^a'^^ + Zi : 008/3$ :^ (4.7) 



^(x) =: cos («$ + 5) : , g = -X. (4.8) 

We recall that 

- = — , m, n G N (4.9) 

p n 

and we denote the soliton mass in the unperturbed theory by Mp, as in Subsection pT] . 
Now it is straightforward to apply the general formulae summarized above, however, some 
care must be taken since in the unperturbed theory (A = 0) the vacuum and the breather 
states come in n-fold degenerate multiplets reflecting the folded nature of the limiting SG 
model. When the perturbation is switched on (A ^ 0) it not only shifts the vacuum energy 
density and the masses but also removes (possibly part of) the degeneracies. 

Note that it is possible to interchange the role of the two cosine terms and therefore 
there are eventually two perturbative regimes of the model. We present the formulas for 
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the above assignment; it is a trivial matter (by suitably shifting the field and redefining 
the S phase parameter) to get the formulae corresponding to the other possibility. We also 
frequently use the parameter 



P 



2 ■ 



(4.10) 



4.2.1 Vacuum energy densities 



First we discuss the shift in the vacuum energy density. From ( p.5|) we see that the model 
^0 has n vacua \k) , k = 0, . . . ,n — 1, characterized by 



{kmx)\k) = jk. 



(4.11) 



In infinite volume, no local operator has nonzero matrix elements between \k) and \k') for 
k 7^ k'; therefore, the perturbative formulas can be applied for each state \k) straightfor- 
wardly. This applies to all excited states above the vacua \k) as well (see e.g. the breathers 
in the next subsection). According to ( |4.3|) the vacuum energy density of the kth vacuum 
is shifted by 



SSk 



X 



Using (|A.4| , |A.25| ) we obtain the final form 



-XGaiP) COS 



27ra 



k + 5 



(4.12) 



(4.13) 



where Qa{f3) = Qa\f3)- We remark that one can extract a similar result from the classical 
potential 



V{^) = -/i cos - A cos («$ + S) . 



At A = 0, the minima are 



(0) 



2nk 



V $ 



-fx 



(4.14) 



(4.15) 



Switching on a small A, the minima are shifted to 



A — — sm 



(27ik^ + 5^ +0{\^) (4.16) 



and the new values of the potential at these minima are 



-/X — A cos 



k + 5 



(4.17) 
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^^classical _ _^ ( + ^ ) , (4.18) 



Therefore we obtain 



which is exactly the classical limit of ( [4.13| ) since it can be easily proven from (|A.4|) that 

01 

^ 1 as /? ^ , - = fixed. (4.19) 

P 

4.2.2 Mass corrections 

We calculate the first order corrections to the masses of the first and the second breathers, 
Ml and M2. In general 

where, as before, Mp is the soliton mass in the unperturbed theory A = 0. 

Eventually, there are n copies of each breather: B). is the rth breather living over the 
vacuum \k) , A; = 0, . . . , n — 1 (c.f. 0). In the theory ^0, their masses are independent of 
k, but the degeneracy is split by the perturbing potential. We denote the corresponding 
mass corrections by oMr . The masses of the first and second breathers can be calculated 
using (|¥^ ) and ([0| , and turn out to be 



M,j V (3 / sin^(f)cos(pT) 

(2cos^(f ) -sm^(p,|))exp{-^ j^^'dtjl^} , (4.22) 

(the constant Af is defined in ( |A.14| )). 

One can also do a semiclassical analysis for the first breather mass. We remark that 
the first breather is identified with the particle created by the fundamental scalar field $ 
and therefore its mass is given by the second derivative of the potential V at its minima. 
Explicitly we obtain 

= ^(3 + cos (^^k + 6^+0 (A^) , (4.23) 

which can be matched with the classical limit /3 — > of (|4.21|) using 



j^^classical ^ ^ ^^ 24) 

and A/" ^ 1 as /5 — > following from ( [A.14| ) . 
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4.2.3 S-matrix corrections 



We can also compute the first order correction to the Bi — Bi S'-matrix, using (|4.6| , |A.9| ). 
The result is 

6Su{^) = I co^[—-k + 5\Su{d) . o.p^x (4.25) 



(4.26) 



Ml V P J ' sin^(f) 

sm^ (p7r|) ^ 1 \ sinhi^ 
sin^ (pvr) cosh('i9) + 1 / sinh^ -d + sin^ (pyr) 



where 



smh -(7 — z sm pvr 

is the sine-Gordon -Bi — Bi S-matrix. It turns out, however, that this correction is too 
small and cannot be measured within the precision of TCSA. 



5 Comparison of FFPT with TCSA 

5.1 Extracting the vacuum energy densities and the mass spec- 
trum from TCSA 



TCSA gives the spectrum of the dimensionless Hamiltonian hrcsA ( |3.14|) as a function 



of the dimensionless volume parameter / ( p.l5| ). In our numerical calculations, we use 
the interpolating mass scale M (p.lO| ) as a unit of mass. However, the FFPT results are 



expressed in terms of the scale M^; we must therefore take care of converting between 
the two conventions. We describe first how we extract the vacuum energy densities and 
the various masses from the TCSA data then we compare these 'experimentally measured' 
quantities and the FFPT predictions on a specific example when a/ j3 = 1/2. 

5.1.1 Vacuum energy density 

Concerning the vacuum energy density, our numerical results show that the ground state 
energy is linear in the volume L 

-2'vac = £va.cL (5-1) 

for a wide range of L. Therefore the (dimensionless) vacuum energy density £^vac/^^ can 
be measured as eyac (^o) /^o by choosing some value Iq in this range (typically between 14 
and 20). We extract the shift in the vacuum energy density using the relation 



S^vac evac (^o) f S- 



(5.2) 



12 



where si'^ is the value measured using TCSA in the unperturbed theory at rj = 0. It's 
value is known analytically as well 

— 77 = — tan — , 5.3 

however, we prefer to use the value extracted numerically because it helps eliminating part 
of the truncation error coming from TCSA. 

5.1.2 Masses 

We know that the energy of a particle a of momentum in finite volume, relative to the 
ground state, has the form 

Ea{L) - E,ac{L) =Ma + (c"^^) (5.4) 

where Ma is the mass of the particle and M is some mass scale characterizing the theory. 
As the interpolating scale ( |3.10| ) agrees with the relevant scales at the two extremal points 
rj = and rj = 1, while in between it interpolates smoothly between them, we can take it 
as the characteristic scale of the DSG model. Therefore the simplest way to measure the 
mass is to extract it as the value of the energy above the ground state at a suitable value of 
the volume in a range (the scaling region) where the energy difference Ea{L) — Eyac{L) is 
approximately constant; (however, this volume cannot be too large since then truncation 
errors would spoil the measurement). Therefore to obtain the mass corrections we use the 
relation 

^ = (6.(/0)-e..e(/0))^-^ (5.5) 

where Iq is a suitable value of the volume parameter, ea(/o) is the relevant eigenvalue of 
( PI ) and Mi°^ is the mass measured at 77 = 0. Once again, the mass ratios 

are known analytically, but we use the values extracted from TCSA in order to eliminate 
part of the numerical errors. 

5.2 A special case: a = (3/2 

We present the comparison between the numerical TCSA data and the analytical first 
order FFPT results for a special ratio of the frequencies. The reason for choosing this 
example is that (1) it is the simplest possible case and (2) this is the one that we treat 
nonperturbatively in detail later on. Let us first discuss the model itself. The action is 

A = j dt j dx Qa^$a^$ + /xcos/3$ + Acos(^^$ + 5^^ (5.7) 
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We denote this model by DSG2(/?, S) (sometimes omitting the parameters in the parenthe- 
ses). The field is defined as an angular variable (|3.5| ), in this case with the period 



$ = $ + —/ , /gZ. (5.8) 

Therefore the theory DSG2~° has two vacua, which are degenerate in infinite volume and 
correspond to the classical minima 

$0 = , *i = ^, (5.9) 

thus this model corresponds to a 2-folded sine-Gordon theory SG(/?,2) [^j. The other 
extremal case DSG2~^ has only one vacuum 

K = -J (5.10) 

and is a 1-folded model SG(/5/2, 1). 

In the vacuum sector of SG(/9, 2), every state essentially appears in two copies, corre- 
sponding to the two diff'erent vacua. In infinite volume, these copies are degenerate, but 
the degeneracy is lifted for finite L by the tunnelling between the two vacua. In particular, 
the vacua and the (neutral) one-particle states are split by an amount decreasing expo- 
nentially with the volume L. Switching on the perturbation removes this degeneracy in an 
interesting way. 

Let us first discuss the case 6 = 0. Switching on A, the two vacua have different 
vacuum energy densities and thus the leading term in the split between the two ground 
states becomes a term linear in L. To first order in A, we obtain from ( |4.13 ) 



(L) - Ei^J,{L) = 2Xg^/,{P)L + ... (5.11) 

where E^^'^\L) are the energies of the two vacuum states (Casimir energies) as functions 
of L and the ellipses denote terms which fall off exponentially with the volume. Since 
the vacuum energy density is universal for all excitations lying above the same ground 
state, the space of states separates into two parts which manifest themselves as lines with 
two different slopes in the large volume limit. In infinite volume, the Hilbert space of the 
system eventually reduces to the ones with the smaller slope, as the other set of states will 
lie at infinitely high energy above them. For this reason we call the states with the larger 
value of the vacuum energy density "runaway" states. 

Of course, higher order terms modify the linear dependence of the slope on A. In 
the other extremal point DSG2~^ one has only one vacuum, but, as the corresponding 
sine-Gordon model is more attractive, it has more breather states. This is where the 
interpolating mass parameter ( |3.10|) comes in handy: as the Hamiltonian (p.l4| ) depends 



smoothly on r], one can run a 'movie' in r/, diagonalizing the Hamiltonian in TCSA for 



different values of rj from to 1. As shown in Figure 5.1, the energy of the set of runaway 
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states first moves upwards, but then turns back and comes down to form the spectrum of 
the endpoint DSCr^ = SG(/?/2, 1). 

We wish to note that the runaway states do not eventually cross the other lines; rather, 
there are line avoidances which is a typical pattern for states that are metastable in finite 
volume and eventually decay ||10| (see also Q for similar results for the off-critical Ising 
model in magnetic field). 

Let us now discuss what happens for general 6. Specializing (|4.13| , |4.21| , |4.22| ), we obtain 





-5^(1) = - 










5m(°) = 




A6^/3/2(/5) 



Ds6 (5.12) 

cos 5 exp < / dt > (5.13) 

^\ Trio sintj ^ ' 

.2 3cos2(f) -1 \ 2 r . t \ . . 

^cos6 exp<^ / dt M5.14 

cospvr TT Jo smt J 

Note that at 5 = | the 0(A) corrections vanish and in particular the two vacua become 
degenerate. We shall show shortly that this is also true nonperturbatively, at least as long 
as rj is smaller than some critical value ?7crit- 

Now we proceed to the comparison between FFPT and TCSA, using plots. The first 
order corrections are linear in the coupling constant of the perturbing operator, so to use 
a convenient and dimensionless parametrization, we define the following combination: 

A = AM^"" . (5.15) 

Figure ^]2| shows the results for the correction to the vacuum energy density 5£^vac/^^ for 
6 = and /3 = iy/n/S, in units of the soliton mass M^. The agreement is spectacular 



for a quite long range of the coupling. Similarly, Figure |5.3| shows the breather mass 
corrections for the four particles Bf'\ Bf\ and at the same values of (3 and 5. 
For comparison we have also plotted the classical prediction for 6Mf'^ . Finally, Figure |5T^ 
presents the dependence of the mass corrections on the phase parameter 5 at a suitably 
small value of A. (Note that here, to ensure better visibility, we changed the scale on the 



vertical axis with respect to Figure As a result, the deviation between the first order 
FFPT and the TCSA results is much more pronounced). 

One can also check the agreement between TCSA and FFPT for a perturbation theory 
in the coupling fj,, i.e. around rj = 1. The results are similar, except that the validity range 
of first order perturbation theory is smaller than around 77 = 0. The dependence of the 
first and second breather masses on r] is shown in Figure [575| . 

To sum up, the comparison between TCSA and FFPT gives a convincing evidence for 
the validity of both approaches in the perturbative (77 ~ or 77 ~ 1) domain. This is of 
utmost importance because for the purposes of further investigations we shall continue to 
use TCSA inside the complete interval < r] < 1, where there are no ways to check its 
validity by any independent method (except some semiclassical analysis and qualitative 
arguments). 
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(b) T] = 0.2 
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(e) V = 0.8 (f) r? = 1 



Figure 5.1: Change of the spectrum as rj varies from to 1 in DSG2 for j3 = 4v^/3 and 
5 — 0. The energies are normalized by subtracting the ground state contribution. The 
plots show the first 20 states (including the ground state, which has constant zero energy 
in this convention). 
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Figure 5.2: The correction to the vacuum energy density as a function of A 
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Figure 5.3: Breather mass corrections as functions of A 
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Figure 5.4: Breather mass corrections as functions of 5 at A ~ 0.01 and (5 — 4v^/3 



Figure 5.5: The dependence of the first and second breather masses on r] running from 
to 1 in DSG2(4\/7r/3, 0). The dots are the numerical TCSA data while the lines are the 
perturbative results obtained by using FFPT around the two endpoints 77 = and rj — 1. 
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6 Phase transition in DSG2: generalities 



We now turn to examining the phase diagram of the theory DSG2 as rj varies from to 
1. We have seen that the theory has two degenerate vacua at 77 = 0, while there is only 
one ground state at 77 = 1. This hints at the possibility of a phase transition happening 
somewhere in between. However, for generic values of the parameters the degeneracy 
between the two vacua is lifted as soon as 77 > thus the end point 77 = is singular in 
this sense. 

However, as we have already pointed out (subsection |5]^) the first order corrections to 
the vacuum energy densities ( |5.12| ) vanish exactly when 6 takes a value on the boundary 
of its fundamental range 



From now on we restrict ourselves to this special value of S and by a slight abuse of notation 
denote this model by DSG2. 

In this section we present an argument that the energy difference between the two 
vacua eventually vanishes to all orders in A and that there is a critical value rj = rjcrit at 
which there is a second order phase transition. Since all perturbative corrections to the 
vacuum energy difference vanish, this transition is entirely in a nonperturbative regime and 
form factor perturbation theory is not applicable. We first present a semiclassical (mean 
field/Landau-Ginzburg) analysis, extending the arguments given by Delfino and Mussardo 

and Fabrizio et al. ||^, and discuss the signatures of first and second order phase 
transitions in finite volume. In the next section we use TCSA to verify the predictions of 
the semiclassical considerations and in particular to establish the second order nature of 
the phase transition and its corresponding universality class. 



6.1 Landau-Ginzburg analysis 

We start our analysis with the classical potential 

V{^) = -/icos/5$ - Acos(a$ + 5) 

= -/icos/3<l> + Asin 1^^$^ . (6.2) 

For definiteness, let us suppose that /x and A are positive without losing generality. As the 
potential is periodic under $ $ + 47r//3, the fundamental range of the field variable can 
be taken as 

— <*<^. (6.3) 
Note that the potential also has the discrete (Z12) symmetry 

T : $ ^ -— - $. (6.4) 
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Figure 6.1: The behaviour of the classical potential V{^) for A < 4/i and A > 4/i 

Then a simple analysis establishes that the potential shows the characteristics of a second 
order phase transition at the point A = 4/i (see Figure ^Tj). Since classically the dimensions 
of /i and A 2, the definition ( |3.11|) of r] reduces to 

/i"^ + A"^ 

and so this point corresponds to the following critical value of r] 

16 

^cnt(/3 = 0) = - = 0.941... . (6.6) 



In view of the Z2 symmetry (|6.4| ) one would like to conclude that this is a second order 
phase transition point in the Ising universality class. However, there is a caveat here, as 
discussed by Fabrizio et al. 0. 

The essence of the problem can be summarized as follows. The potential V{^) is 
renormalized at the quantum level. To determine the phase structure, one needs to look 
at the effective potential. Due to the periodicity and the Z2 symmetry, it can be readily 
checked that the only possible new terms are of the form 

cos (r/3$) , r e Z or sin (r/3$) , r G Z + ^ . (6.7) 
The conformal dimensions of these operators are 

A? = ^ . (6.8) 

and so as one increases r they are less and less relevant. The first two possibilities are 
particularly interesting. Fabrizio et al. ||2| claim that while sin (^$) has no effect on the 
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nature of the phase transition (apart from modifying the position of the critical point), the 
term cos (2/3$) which is a relevant perturbing operator for < 2-n can (depending on its 
coefficient) turn the second order phase transition into a first order one. 

Let us analyze the effect of such operators now. We take as our quantum effective 
potential the expression 

V;ff($) = -/i cos /3$ + A sin + sin ^<l> + v cos 2(3^ . (6.9) 

K and v can eventually be expressed in terms of A and /i, but their precise expressions are 
unknown. The conditions for the existence of a second-order phase transition point are 

,>0, (6.10) 

since the Z2 symmetric extremum is at $ = — ^. 
Explicitly, we have the conditions 

-4/i + A - 9k - 16z/ = , 16/i- A + 81k + 256z/ > . (6.11) 

Putting K = z/ = we recover that the phase transition point in the classical potential is 
at A = 4/i and the second derivative is positive (we assumed that A and /x are positive) . 

Suppose now that k and v are independent parameters and that we vary only A. The 
critical point can be read off from the first condition in (|6.11| ) : 

Ac = 4/i + 9ft + 16z/ . (6.12) 

Therefore the condition for a second order transition is 



, v^^m 



/i + 6K + 20z/ > . (6.13) 

We can draw three conclusions. First, there is ample space for second order phase transition 
even in the presence of the higher frequency terms in the effective potential. Second, the 
term sin^$ is also capable of inducing a first order transition, since a large negative k 
is enough. Third, it is also possible that higher order derivatives vanish at the transition 
as well, giving tricritical and in fact even tetracritical points (even higher frequency terms 
open the possibility to multicritical points of any order). These observations can also be 
confirmed by simply plotting the effective potential for different values of the parameters. 

However, for the double sine-Gordon model (|2.1J ) the coefficients of the higher frequency 
terms (e.g. /t, v) are fixed by the dynamics. Therefore this model lives on a specific surface 
in the multiparameter space of possible quantum effective potentials that respect the Z2 
symmetry $ — > — ^ — $ and the periodicity $—><!> + 47r//3. As we have already shown, 
the phase transition point is out of the reach of perturbative techniques and therefore one 
needs to recourse to some nonperturbative analysis to solve the problem. In the special 
case when p = 1 (/5^ = 47r), by using a mapping onto a generalized Ashkin Teller model, 
Fabrizio et al. were able to show the existence of a second order phase transition. In 
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the general case, when < p < 1, since the model is nonintegrable, the only available 
candidate is TCSA (the upper limit is to avoid any UV problems with TCSA). 

Finally we present the characteristic picture of a first order transition (see Figure |6^ ). 
Since we disagree with the claim made by Fabrizio et al. that sin cannot induce a 
first order transition, we present a case in which it does. We fix the following values of the 
parameters 

2 

K = --ii , z/ = 0, (6.14) 



and let A vary from upwards. At A = 0, we have only the two vacua in which Z2 is spon- 
taneously broken. At A = |/i the symmetry preserving vacuum appears as a metastable 
state. For A = y/i it becomes degenerate with the symmetry breaking vacua: this is 
the so-called coexistence point where the three vacuum states exist simultaneously in the 
theory. Then for A > ^/x the symmetric ground state becomes the true vacuum and the 
two other are metastable ones which disappear at A = 

6.2 Signatures of 1st and 2nd order phase transitions in finite 
volume 

In order to establish the nature of the phase transition we shall examine the spectrum of 
the model in a finite volume. To prepare the ground we now discuss the signatures of both 
types of phase transitions in the finite volume spectrum. 

6.2.1 2nd order phase transitions 

For definiteness we specify the symmetry of the transition to be Z12, although the discussion 
can be easily generalized. In the broken phase, there are two degenerate vacua in infinite 
volume. For finite L the degeneracy is lifted by tunneling effects through the barrier 
separating them and the lowest lying becomes the ground state, while the other one is 
the first excited state. The split between the ground states vanishes exponentially as 
L — > 00 because the height of the barrier is proportional to L. For L — the two ground 
states tend to some eigenstates of the conformal Hamiltonian and therefore the split is 
proportional to . 

In the unbroken phase there is only one vacuum. The first excited state is generally 
some massive one-particle state and has a finite gap M. over the vacuum. Therefore the 
difference between the ground state and the first excited state tends to as L — cxd in 
the same way as indicated in eqn. (|5.4| ). 

Let us call the parameter in which the transition happens 77crit, in analogy with the 
DSG model. Then the characteristic behaviour can be summarized if we look at the mass 
gap as a function of t] and suppose the unbroken phase is 17 > r/cHt while the broken 
one is ?7 < r/crit- We have the following behaviour 

M{ri) = , r]< r/crit 

Miv) > , V>Vcrif (6.15) 
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Figure 6.2: Illustrating a first order transition with k — —2/5/1 and u — 
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At exactly rj = rjcrit by changing L from to oo we interpolate between the UV c = 1 CFT 
and an IR fixed point one, which consists of an IR CFT and a massive one. Therefore 
certain states in the spectrum have the asymptotic (large L) behaviour 

E*(L)-Eo(L) = ^^^i^5L±M + o(L-i-) , e>0, (6.16) 

where the subscript denotes the ground state (the lowest lying state in finite volume). 
These states fall into conformal families of the infrared fixed point characterizing the tran- 
sition with conformal dimensions Aj^. As L — > oo they decouple from the rest of the 
spectrum which remains massive. The identification of such states allows the extraction of 
the characteristics of the infrared fixed point CFT as we shall see shortly. 

6.2.2 1st order phase transition 

A first order phase transition is always characterized by the existence of metastable states. 
As we have seen in Subsection |5]^, the existence of such states in finite volume implies 
that the space of states separates into excitations over the true ground state and those 
over the metastable ones which we called "runaway" states. The two sets of states have 
in general different bulk terms which coincide only when one tunes the coupling rj to the 
coexistence point ?7coex- At exactly this point the three vacua have identical energy densities 
and passing through this value of rj the symmetry breaking pair of vacua interchanges with 
the symmetric one. This is a very characteristic behaviour which one would think is very 
easy to distinguish in the finite volume spectrum. 

However, first order phase transitions have a characteristic strength and as is well- 
known, weak first order phase transitions are hard to distinguish from second order ones. 
To see this more concretely, let us denote the three vacua by $o for the symmetric one and 
$-1- for the symmetry breaking ones, respectively. There are two other interesting points in 
rj: the lower critical value ?7icrit where $o appears, and the higher critical value r/hcrit where 
$-1- disappear as we increase rj. In the case of a second order phase transition 

^Icrit = ??hcrit • (6.17) 

If the difference r/hcrit ~'7icrit is small, then when tuning rj numerically (or in an experiment) 
one may miss the existence of the metastable vacua. Similarly, one can look at the difference 
between the vacuum energy densities V {^o) — V ($±) at the lower and upper critical points. 
If it is small, then to detect the fact that the space of states separated into two sets of 
states with different linear terms in their energies one needs to go to very large values of 
the volume where truncation errors in TCSA become large. Therefore it is clear that one 
can never eventually distinguish a sufficiently weak first order transition from a second 
order one, neither numerically nor experimentally. An exact solution of the model would 
of course disentangle the problem, but given that the theory is not integrable this is too 
much to hope for. 
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7 The phase diagram in the case | = |, ^ = \ 

In this section we use TCSA to verify the classical predictions for the phase transition, and 
to clarify the order and the universality class of the transition. In order to check for the 
transition we need TCSA data over a wide range of the volume parameter /. Therefore, 
in order for TCSA to be useful, we have to restrict our investigations to the domain 
(3"^ < 47r where the method does not have UV divergences (see the discussion at the end 
of Subsection |3^) . It turned out that even in this range we cannot attain the necessary 
precision everywhere and to get a convergent enough numerical method we had to restrict 

< Svr/S. For the 77 = sine-Gordon theory, this corresponds to a regime with two or 
more breathers in the spectrum. 

To examine the nature of the phase transition we ran the t] 'movie' with 5 = tc/2 at 
various values of /?. The results show no trace of a first order phase transition - a case 
is depicted on Figure |7]T[ Indeed, for the case shown (/? = 8\/7r/7) the spectra can be 
classified into three classes: (1) for rj < 0.8 one can clearly see the doubly degenerate 
ground states and the also degenerate first massive states above them; (2) for 0.9 < r] < 1 
the spectrum is massive, but there are no degeneracies among the lowest lying states and 
(3) in the transition domain, for rj ~ 0.84, the structure of the spectrum changes, there 
are no lines tending to a constant, instead they decrease in the whole volume range and 
the degeneracies of the ground states and the first excited states are removed. Since at no 
place could we see the linearly rising set of runaway states that would correspond to the 
presence of metastable vacua with a different bulk energy, we analyze the data by looking 
for a 2nd order phase transition at some r] = r/crit- (See however the remarks at the end of 



Subsection 3.2.2). 



7.1 How do we find 77 = r/crit? 

We determined the critical value of r] in the following way: we tuned rj in the transition 
region by looking for whether the €2(1) — evac(0 difference of the second excited state and 
the ground state continues to decrease along the (in between enlarged) complete / range. 
The meaning of this criterion is clear: this difference describes both at rj = and at 
T] = 1 the mass of some low-lying breather of the limiting SG theories; thus its continuous 
decrease at a specific rjcnt implies that this breather became massless. Since this happens 
only at this particular value of rj, we think that this criterion is better than the one based 
on the appearance of a gap between the two ground states (as the gap is there for all 
V > ^crit)- Nevertheless, for consistency, at 77crit chosen, ei(/) — evac(0 should decay slower 
than exponentially with I; in fact, recalling eqn. (|6.16| ), it should be 0(1//) for large /. 



The result of this search is shown on Figures |7^ - |7^ for three different values of p. 
All three plots were taken at the 'best' rj values according to our criterion. To separate 
more clearly the various ei{l) — eo(/) as functions of / on all three plots we show these 
differences multiplied by //(27r). (This also makes them more suitable for the subsequent 
analysis). Here eo(0 denotes the lowest eigenvalue (vacuum energy) and i runs from to 
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Figure 7.1: Change of the spectrum as rj varies from to 1 in DSG2 for /? = 8\/7r/7 and 
S = n/2. The energies are normalized by subtracting the ground state contribution. The 
plots show the first 8 states. 
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Figure 7.3: [ei(l) - eo(l)] l/(27r) at /3 = A^/7^/3 and 77 = 0.850 



27 



1 r"T""r"T""r"-f""r"r"T""T""r 

20 40 60 80 100 120 

Figure 7.4: [e,(/) - eo(/)] //(27r) at p = 8^/5 and r] = 0.838 

7. The horizontal lines on the plots correspond to the possible values of A^^ + Ajj^ if we 
assume that the IR CFT is the c = 1/2 Ising model (see eqn. (|6.16| )), when the spectrum 
of Ajp^ + Aj^ is determined by the Ising primary fields and their descendants. 

In all three cases we see that the first two lines - that correspond to the first two excited 
states above the ground state - qualitatively fit rather well to the Ising model's prediction. 
For the higher excited states the agreement is not so spectacular but even for them the 
sequence of lines and the absence of degeneracies (which on the Ising side follow from the 
null vector pattern) match that of the Ising prediction. For the two smaller /3 values it is 
clear that also some of the higher excited states would reach their corresponding plateaus 
by going to even higher volume^. 

To make this agreement more quantitative we have to describe two types of effects. On 
the one hand we need the spectrum of states on the critical trajectory (i.e. at rj = r/crit) 
in a finite volume more precisely than in eqn. ( |6.16| ), and on the other we need a clear, 
quantitative picture of the deviations from criticality since our best rj may not equal ?7crit 
(and TCSA errors may also be involved). 



^These plots show that the higher the value of (3 the larger are the TCSA errors. The spectra were 
computed by using a conformal energy cut Ecut ^15—16 which resulted in 4.5 — 5.5 x 10'^ states in the 
UV conform Hilbert space. 
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7.2 Theoretical predictions for the form of the spectra 

We describe these effects by assuming that the IR CFT is the critical Ising model, i.e. by 
using the Hilbert space and operator content of the Ising model. 

Let us concentrate first on the finite volume corrections at rj = rjcrit- In this case when 
changing I from zero to infinity we interpolate from the UV conformal fixed point to an IR 
one, at least if we neglect the coupling between the massive and massless modes in the IR. 
The massive and massless modes decouple exactly at Z = oo and we shall call the subspace 
of states asymptotically corresponding to the massless ones the scaling sector. Since one 
can flow into a fixed point only from an irrelevant direction, in this approximation the 
interpolating theory may be considered as an irrelevant perturbation of the conformal 
Ising model with (the non renormalizable) action 

A — Aising + g J (i^-2(^(-2, ^) + higher dimensional fields) . (7.1) 

The various directions correspond to various choices oiip{z, z), which in turn are expressed 
in terms of the irrelevant (spinless and non vanishing) descendants of the primary fields. 
Since the DSG2(/9, 7r/2) model exhibits the Z2 symmetry for all values of r], ijj{z,z) may 
contain only the descendants of the energy operator dde{z, ^)+higher dimensional fields and 
of the identity operator TT+higher dimensional fields (where T{T) is the (anti)holomorphic 
component of the stress energy tensor of the Ising model) . Recalling the generic form of 
the TCSA Hamiltonian, 

we see that for large but finite / the leading corrections to E^{L) — Eq{L) come from the 
least irrelevant operator: 

9'jr 

e^il) - eoil) = y (A+ + AfJ + A^l'-'"^^ + ... (7.3) 

(Here the constant A may depend on the state in question). Thus, in our case, the leading 
correction takes the form 

6^(0 - 6o(0 = -(A+, + A,!,) + A^r' + ... (7.4) 

if li' is given by the derivative of the energy operator. (If this term were absent then TT 
would give a correction A,!,/^'^). 

Numerically we can probably never tune r] exactly to rjc, and even if we could the 
TCSA errors - coming from truncating the Hilbert space - would drive away from the 
critical trajectory connecting the UV and IR conformal theories. From the side of the 
Ising model this means that some relevant perturbation is also switched on, and we should 
also take its effects into account when discussing the large / spectrum. In the Ising model 
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there is only one relevant perturbation compatible with the Z2 symmetry, namely when the 
perturbing operator is e{z, z). The presence of a relevant perturbation 7e(l, 1) in the TCSA 
Hamiltonian leads to a correction B^, or B^, + C^,l to e^,(/) — eo(OQ- Thus, summarizing, 
we expect, that in a large but finite volume range the energy spectrum can be described 
by the following: 

eS) - eo{l) = B, + —A + Ar' + + (7.5) 

In the following we fit the TCSA data using this expression. We keep all these terms since 
we are at an intermediate volume range, where all of them can be equally important. A 
little bit more abstractly we can say that on the "best rj" trajectory we come into the vicinity 
of the IR fixed point, where the features of this CFT determine the energy spectrum, but 
eventually, because of the relevant perturbation, we are driven away and fiow into a massive 
theory. Nevertheless, if in the intermediate volume range these expressions describe the 
spectra well, - in particular if the measured values of Bi (or Bi and Ci) are small and that 
of Di are similar to the ones predicted by the Ising model - then this confirms the presence 
of a 2nd order phase transition in the Ising universality class. We call this intermediate 



volume range the scaling region which may (and in fact does, see Figures |7^, |7.3| and [L 
change from state to state. 



Let us pause for a short discussion of the plateaus apparent in Figures |7^ , [0| and 
Note that if we had the exact finite volume spectrum (i.e. without truncation errors) then 
the scaling functions 

f (e.(0-eo(0) (7.6) 

ZTT 

would tend to a constant as / 00 for 77 = r^crit. In this sense the IR fixed point lies 
infinitely far away and the above functions approach asymptotic "plateaus" at / = 00. 

In the real TCSA situation, the scaling region corresponds to the plateaus in the nu- 
merical scaling functions. Tuning rj closer and closer to the critical value these plateaus 
become longer and fiatter and in addition they move to higher and higher values of the 
volume /. 

Experience with TCSA for perturbed c = 1 CFTs shows that truncation errors generally 
decrease if one decreases /? (or equivalently, increases the compactification radius of the 
conformal boson field). As a result of the above considerations, we expect that the plateaus 
are longer and fiatter for smaller /3 and that they start at higher values of I. This is in fact 
what we observed in practice (see Figures ^]3| and [7^ ). 

One can also make a theoretical prediction for the form of the energy eo(0 of the ground 
state as well. Using well-known facts about the off-critical Ising model ||ll| and arguments 
similar to the above one expects the following behaviour in the scaling region 

eo(0 = -TT7 + oo^' + bo + Col + Co/ log / + . . . (7.7) 
bl 



^The second term seems to be dominant over the first one, however its coefficient - which in a pCFT 
framework is of second order - in general is much smaller than the constant, thus for a large but finite 
range of ^-s they may be equally important. 
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State 








a 


i = 1 


0.125 ±0.002 


-1.2 ±0.3 


-0.0034 ±0.0001 


1.45 X 10^^ ± 6 X 10"'^ 


i = 2 


1.04 ±0.03 


-128 ±8 


0.0017 ±0.0018 


9 X 10-6 ± 5 X 10-6 



(a) P = 80F/7, T] = 0.866 



State 




A. 


B, 


a 


i = 1 


0.1312 ±0.0009 


-1.35 ±0.03 


-0.0061 ±0.0003 


2.6 X 10-5 ±6 X 10-6 


i = 2 


1.03 ±0.02 


-74 ±3 


0.005 ±0.002 


1.4 X 10-5 ±9 X 10-6 


(b) /? = 40F/3, 77 = 0.850 


State 




A. 


B, 


Q 


i = 1 


0.1455 ±0.0005 


-0.89 ±0.05 


-0.0145 ±0.0005 


2.24 X 10-^ ±9 X 10-6 


i = 2 


1.09 ±0.02 


-40 ±2 


-0.0015 ±0.0031 


1.1 X 10-^ ±2 X 10-5 



(c) P = 80F/5, ?7 = 0.838 



Table 7.1: The results of fitting ( |7.5|) to the first two excited states for various values of (3 
at the estimated critical value of r] 



where c = 1/2 is the central charge of the Ising model and the logarithmic term appears 
as a result of a "resonance" between the linear (bulk) vacuum energy contribution and the 
perturbative corrections coming from the relevant operator e. The term CqI log I is universal 
for all the excited states in the scaling sector and so disappears from energy differences 



like (7.5). 



7.3 Numerical results 

7.3.1 Scaling dimensions and UV-IR operator correspondence 

One can fit the function ( [7. 5]) to the TCSA data in the scaling regions plotted in Figures 



[7^ , fTdj and [7^ . The results of the fits are shown in table [7T . 

From the conformal data of the Ising fixed point we know that the first excited state 
should correspond to the spin operator a with conformal dimensions = 1/16, while 
the second to the energy operator e with conformal dimensions = 1/2. Therefore we 
expect 

Di = ^ = 0.125 , 1^2 = 1. (7.8) 
8 

The data presented fit quite well with these predictions. Note that the agreement becomes 
better for smaller /? where the TCSA errors are smaller. We mention that the errors 
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presented in the tables come from the fit procedure and do not contain the truncation 
errors which are generally much larger. 

Another piece of information is the UV-IR operator correspondence. In TCSA we know 
the operator in the UV theory which creates the excited states. By identifying these states 
in the language of the IR theory, one can establish a correspondence between the UV and 
the IR operator content. In particular, the first and the second excited states are created 
by the operators 

, . 

cos-$ and sin (7.9) 

in the UV, respectively. In the IR these states fiow to ones corresponding to the Ising 
operators a and e. This is in full accordance with their parity under the Z2 symmetry ( |6^ ) 
and the fact that a is an odd, while e is an even operator in the Ising model. In addition, 
both the operator sin |$ in the DSG2 model and the operator e in the Ising case correspond 
to the "temperature" parameter deforming the model away from criticality. Furthermore, 
this conclusion also agrees with the conjecture made by Fabrizio et al. ||2|. 

One could also try to extract the central charge c using the predicted form of the 
vacuum energy function ( |7.71 ). It turns out, however, that the fit is not sensitive enough 
to the value of c — in a least-squares algorithm the main contribution to comes from 
the linear and / log / terms and c can only be determined with 30 — 40% accuracy, within 
which we find agreement with the Ising value c = 1/2. 

7.3.2 Analysis of the phase diagram 

Concerning the phase diagram, our first observation is that all evidence points to a second- 
order phase transition for all values of /3, which happens at some critical value r7crit(/5) 
dependent on (3. As we already discussed in Subsection |6]^, there is always a possibility 
of a sufficiently weak first order transition. The only statement one can make is that the 
phase transition is second order up to the precision within which the coefficients Bi and 
Ci in (|775|) vanish at the (approximate value of the) critical coupling. Given the data in 
table ^T] and the fact that the normal precision in TCSA is of the order 10^^ — 10^^ for 
the values of the volumes used (this can be estimated e.g. calculating an integrable case 
when the exact values of the excited state energies are known) we can state that the phase 
transition is second order up to the precision of the TCSA method. 

As we have already seen the second order nature of the phase transition is controlled 
by inequalities like ( |6.13| ). Since a very weakly first order transition requires a very weak 
breaking of these inequalities this in turn strongly constrains the parameters in the effective 
potential. However, the coefficients in the effective potential are uniquely determined in 
terms of A and /i and therefore such a fine tuning is highly unlikely to happen. 

It is expected that for 327r/9 < < Svr the transition is always second order since 
then all the possible corrections ( |6.7D to the effective potential correspond to irrelevant 
operators. Therefore it is only necessary to consider the case < < 327r/9. As TCSA 
is convergent for < 47r it can be used to examine this range of couplings and all the 
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values of (3 were from this range. Unfortunately, TCSA does not converge very well for 
> 87r/3, so there is a gap in the values of up to 327r/9 which is not accessible. 
One also expects a second order transition around /? = since in this case the theory is 

semiclassical and the contributions to the effective potential are expected to be very small. 

We also know the exact value 



r/crit(/3 = 0) = — = 0.941... . (7.10) 

Therefore we conjecture that the phase transition in DS(jI is second- order for all values 
of (3 in the exact theory too. 

One can prepare a phase diagram in the plane. For 77 = the theory is in the 

broken symmetry phase, while for 77 = 1 the Z12 symmetry is restored. The phase transition 
line is described by a function r] = ?7crit(/3). One can also calculate the value of this function 
at the point (3 = \/8tt. For l3 > a/Stt the term cos (3^ becomes irrelevant and the theory is 
always in a symmetric phase. Therefore we expect that 

Vcnt{(3 = V8^) = 0. (7.11) 

Since the transformation /3 —>■ —[3 can be compensated by a field redefinition $ — $, 
one also has 



Vcrit {-(3) = Vcrit iP) 



(7.12) 



Figure |7^ shows the phase diagram of the theory. The numerical data are the TCSA 
results for ?7crit at particular values of /? (in addition to the ones in Table fTTI , there is also a 
point with f3 = 4-y/7r/7 where the measured critical coupling is ?7crit 0.916 ±0.002). The 
continuous line is an (even) polynomial interpolating the numerical data together with the 
two known values at /3 = and v^Stt, which is shown in order to get an idea of the phase 
transition line. Below the transition line the theory is in the broken phase, while above 
the line the Z12 symmetry is restored. 



8 Conclusions 

In the present paper we examined the two-frequency sine-Gordon model, using a combina- 
tion of various (analytic and numerical) approaches. Using form factor perturbation theory 
(FFPT) on the one hand and truncated conformal space approach (TCSA) on the other, 
we first verified the consistency of the two by comparing the results in the two perturbative 
regimes. Based on the agreement between the two methods, we continued to use TCSA in 
the nonperturbative regime. 

Then we set up a framework for examining phase transition using the finite volume 
spectrum extracted from TCSA. We discussed how to distinguish between first and sec- 
ond order transitions, and in particular how to extract the characteristic quantities of an 
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Figure 7.5: The phase diagram 



eventual IR fixed point. To our knowledge this is the first case when TCSA was used for 
locating and characterizing a nontrivial fixed point. We would like to mention that TCSA 
was used before in the similar context of integrable massless flows in order to check results 
from Thermodynamics Bethe Ansatz (TBA) close to the UV CFT (small values of the vol- 
ume /) P|. Here, on the contrary, the TCSA method was used for finding the critical point 
at which the flow becomes massless and then analyzing the intermediate volume range in 
which the theory is governed by its IR flxed point. In addition, in our case there are no 
alternative nonperturbative methods since the theory is nonintegrable. 

Using this framework, we examined the prediction by Delflno and Mussardo of a 
second order phase transition, in a particular case when the ratio of the two frequencies 
was 1/2. Contrary to a possibility raised by Fabrizio et al. ||2| we argued that this phase 
transition is second order for all possible values of the sine-Gordon-like parameter /?, at 
least within the precision attainable by TCSA. We examined the spectrum in detail and 
showed that the lowest lying states match the pattern of an Ising type fixed point, both 
qualitatively and quantitatively. We also examined the UV-IR operator correspondence 
and found an agreement with the conjectures made by Fabrizio et al. 

In addition, by a simple mean field analysis we have shown that contrary to the claims 
made by Fabrizio et al., all lowest order correction terms to the effective potential are able 
to induce a first order transition, although the numerical studies suggest their coefficients 
are such that they do not affect the second-order nature of the transition in this particular 
case. 

In this paper we restricted ourselves to the investigation of the vacuum sector. It could 
be interesting to see the effect of the phase transition in the topologically nontrivial sectors 
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as well. Another open question is to examine other possible frequency ratios and to extract 
the phase diagram for the general case. Yet another direction for further study is the case of 
irrational frequency ratios: as we already discussed in the paper, the perturbative approach 
(FFPT) is still applicable, but there are no nonperturbative alternatives at present. Finally, 
as the model has potentially interesting condensed matter applications, it is also a challenge 
to use the above framework to extract some quantities of interest in that context. 
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A Vacuum expectation values and breather form factors 
in (folded) sine-Gordon theory 

Here we list the vacuum expectation values and form factors necessary for performing the 



form factor perturbation theory calculations in the main text, using the results in ||T2|, |T3 
We define our conventions for sine-Gordon theory with the action 



AsG = j dt j dx Qa^*^'^^ + /i : 003/3$ :^ 



The vacuum expectation value of the exponential operator 

Va{x) = : exp (m$(x)) : 

= (o|K(o)|o> 



(A.l) 



(A.2) 
(A.3) 



is of the form 12 



4tt 



87r-/32 



X exp 




2 sinh (^gt) sinh (t) cosh (^1 - g) t 



477 



(A.4) 
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provided we normalize our operators in the following way: 

1 



[0\Vaix)V^aix')\0) 



'T'' 47r 



as |x — x'l — > 



(A.5) 



\x — X' 



Here M denotes the soliton mass and |0) is the vacuum in which the sine-Gordon field $ 
has vanishing vacuum expectation value. 

Let us use the following notations for the breather form factors 



where Bk denotes the kth breather with mass 



7T Aj 

Mk = 2M sin ( —p 



P 



Then we have from 



and 



(A.6) 



(A.7) 



(A.8) 



^<k<j<4: 



X 



[a]^ + {ala^ + al) JJ {xi + Xj) ' I 

l<i<j<4 J 



where we used the following notations 



Xk = exp 



sm (pvr^ 



- pn I P7C 

A = 2 cos — \ / sm — exp 
2 V 2 ^ 



sm pn 

■P" dt t 







27r sin t 



44 4 
(Ji = ^ , 0-4 = J]^ Xi , 0-3 = 0-4^ x''^ 



i=l 



i=l 



i=l 



(A.9) 



(A.IO) 
(A.ll) 

(A.12) 



and in the strip — 27r + pvr < Qmi} < —pn (note that in the attractive regime p < 1) the 
function R{{}) is defined as 



„ . , , sinh(t) sinh(»t) sinh((l + p)t) ,0 
R{^) = ATexp <{ 8 / — — ^ sinh^ 



(A.13) 
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with 

M = exp |4 j d^-^^m^^^^!^^^^^±M^ . (A.14) 
R{'d) satisfies the useful identity 

R{^)R{'d±ni) = . , " (A.15) 

smn tf =p sm pre 

which can be used to continue it analytically out of the validity range of its integral 
representation. 

For the mass correction of the second breather we need the explicit form of the form 
factor F22^ . This is not available in the literature, but can be computed from F^^l^ using 
the bootstrap equation 

i Res Fa,...a,Ai (^1, • • • , ^n, + iU^ ' h ^"+1 ' + f ) = (A.16) 

Ti,Fa,...a„2i^l,...,^n,^n+l) • (A.17) 

The fusion angle and the three-point coupling are defined from the B2 pole U^i in 
the Bi — El S'-matrix: 



i{Tl^f g .jj2 1 o /Q^ sinh?? + isinpvr 

n rr2 > ^^-^t^^, ,where Su{^) = —— 

v — lUii smh v — I sm pvr 



and from the relations 

Ul2 = Ti-Ul2 , U^, + 2Ul2 = 2n. (A. 19) 

This way we find 

U^^ = 2UI2 = pTi , (r^^)^ = 2tanp7r . (A.20) 
The residue of R can be calculated from (|A.15| ) 



i 1 

R {-ipn + e) = — tanpTr— — -— (A. 21) 



and finally the B2 — B2 form factor is 



Fif{A,A) = (r?i) 'K22{^)R{^fR{-^-tp7r)R{-^ + tp7c) 



\2 ' 



R{-m{l +p)y 
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Finally we recall that in the perturbative investigation of the DSG model we need the 
vacuum expectation values and form factors in th sectors over the vacua \k) characterized 
by the property 

{k\^{x)\k) = '^k. (A.23) 

Using the following symmetries of the folded sine-Gordon model 

27r 

$^$ + — , (A.24) 

it can be easily shown that 

6?f (/?) = {k\VM\k) = Qa{(i)e'i^' . (A.25) 

The form factors corresponding to states above these vacua can be obtained by substituting 
GaiP) with gi''\p) in the form factor formulae (|A]|, |A]|, |A.22| ) above. 
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